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ABSTRACT 

We present new and archival Chandra snapshot (w 10 ks each) observations of 15 optically 
identified (from the Sloan Digital Sky Survey, SDSS) Type 2 quasai's at z=0.40-0.73. When 
combined with existing X-ray data, this work provides complete X-ray coverage for all 25 
radio-quiet Type 2 quasars with logLpm] > 9.28 Lpj from Zakamska et al. (2003). Two targets 
out of 15 were not detected by Chandra and most of the remaining sources are X-ray weak, 
with nine having less than 10 counts in the 0.5-8 keV band. Low-to-moderate quality spectral 
analysis was limited to three sources, whose properties are consistent with the presence of 
column densities in the range cm ^ in the source rest frame. If the [O III] 

luminosity is a reliable proxy for the intrinsic X-ray luminosity, the current X-ray data indi- 
cate that Compton-thick quasars may hide among « 65 per cent of the SDSS Type 2 quasar 
population (ix,meas/i'x,[om] <0-01); however, since the Type 2 quasar sample is selected on 
[O III] luminosity, the estimated Compton-thick quasar fraction may be overestimated. Using 
archival Spitzer observations, we find that « 50 per cent of SDSS Type 2 quasars appear to 
be obscured by Compton-thick material based on both the ix,meas/i'X,mid-iR (where mid-IR 
corresponds to rest-frame 12.3|a.m) and Lx.meas/i-x.loiii] ratios. We use this information to pro- 
vide an estimate of the Compton-thick quasar number density at z~ 0.3-0.8, which we find 
is in broad agreement with the expectations from X-ray background models. 
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1 INTRODUCTION 

Over the last decade, the quest for the identification of Type 2 
quasarj^ has received renewed interest, and many sources or 
promising candidates of this class have been detected, mainly 
through X-ray surveys (e.g., Norman et al. 2002; Mainieri et al. 
2002, 2007; Fiore et al. 2003; Gandhi et al. 2004, 2006; Ma- 
teos et al. 2005; Severgnini et al. 2006; Delia Ceca et al. 2008; 
Krumpe et al. 2008). Predicted by unification schemes of AGN 
(e.g., Antonucci 1993), they are thought to play an important role 
in many synthesis models that seek to explain the origin of the 
X-ray background (XRB; e.g., Gilli, Comastri & Hasinger 2007). 
On the other hand, the availability of deep radio and Spitzer ob- 
servations up to 24|i.m has allowed to define complementary meth- 
ods to search for and characterize candidate Type 2 quasars (e.g.. 
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obscured (column density Na > 10^- cm^^) Active Galactic Nuclei 
- in X-rays, often characterized by high-ionization, narrow emission lines 
and the lack of broad emission lines in the rest-frame optical/ultraviolet 
spectra. 



Martinez-Sansigre et al. 2005, 2008; Houck et al. 2005; Higdon 
et al. 2005, 2008; Weedman et al. 2006a,b,c; Dey et al. 2008; see 
also Donley et al. 2008 for a review), and examples of the most 
heavily obscured, Compton-thick AGN (with column density above 



n,24 



2005, 2008; Martinez-Sansigre et al. 2007; Daddi et al. 2007; Fiore 
et al. 2008, 2009, hereafter F09; Polletta et al. 2008; Lanzuisi et al. 
2009). 

An alternative approach to X-ray and mid-infrared (mid-IR) 
surveys for the quest of Type 2 quasars consists in following-up, 
in the X-ray band, promising Type 2 candidates selected through 
optical emission lines. In this regard, a careful selection on the ba- 
sis of the [O III]5007A luminosity allowed Zakamska et al. (2003; 
hereafter Z03) to define a sample of 291 Sloan Digital Sky Survey 
(SDSS; York et al. 2000) Type 2 quasar candidates in the redshift 
range ^ 0.3-0.8 (see also Zakamska et al. 2004, 2005, 2006 and 
2008, and Liu et al. 2009 for observations and properties of these 
objects at other wavelengths). These sources are classified on the 
basis of high-excitation, narrow emission lines, without underly- 
ing broad components, and with line ratios characteristic of non- 
stellar ionizing radiation. Standard AGN photo-ionization models 
are able to successfully reproduced most of the emission-line ratio 
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diagrams, at least for the most [O IIl]-luminous AGN in the Z03 
sample (Villar-Martm et al. 2008). 

In this context, the X-ray band is crucial to confirm the AGN 
nature of these sources and assess the presence of obscuration, 
given all the uncertainties related to obscured AGN selection on the 
basis of the [O III] emission line (see ^for extended discussion). 
In Vignali, Alexander and Comastri (2004a, hereafter V04; see 
also Vignali, Alexander & Comastri 2004b), we presented the ba- 
sic X-ray properties of a sub-sample of these Type 2 quasar candi- 
dates, mostly based on ROSAT observations. In Vignali, Alexander 
& Comastri 2006 (hereafter V06), we placed significant constraints 
on the fraction of heavily obscured, possibly Compton-thick AGN 
among the population of SDSS Type 2 AGN, being of the order of 
~ 50 per cent (see also Ptak et al. 2006 for similar findings from 
an X-ray perspective). We note that this result was possible by as- 
suming the correlation between the [O III] and the 2-10 keV flux 
found for Seyfert 2 galaxies (Mulchaey et al. 1994; hereafter M94), 
which is characterized by a significant scatter. Unlike the Type 2 
quasar population found in moderately deep and ultra-deep X-ray 
surveys, at the optical magnitude limit of the SDSS, Type 2 quasars 
are generally easy to study at both optical and X-ray wavelengths, 
as confirmed by the results obtained over the last four years through 
a snapshot strategy with both Chandra and XMM-Newton. 

In this paper we extend the work shown by V04 and V06 by 
presenting Chandra observations for 12 additional SDSS Type 2 
quasar candidates, thus providing a complete X-ray coverage for all 
25 radio-quiet Type 2 quasars with logLjom] > 9.28 L0 from Z03. 
We have also included in our analyses three pointed Chandra ob- 
servations retrieved from the public archive. Coupled with previous 
X-ray observations of AGN drawn from the original Z03 sample, 
the current targets bring the total number of sources with available 
Chandra and/or XMM-Newton constraints to 31 (74 per cent of 
X-ray detections, but see also Lamastra et al. 2009 for recent up- 
dates). 

The work by Z03 has been recently expanded and updated 
by Reyes et al. (2008; hereafter R08), where the [O III]5007A line 
luminosity is computed using both a Gaussian feature and a non- 
parametric line fitting method (see their Section 2.4). More than 
90 per cent of the Type 2 quasar candidates reported by Z03 were 
recovered from this more complex analysis. Since our project was 
originally thought to provide a reliable and complete X-ray char- 
acterization of all the sources possibly classified as Type 2 quasars 
from the Z03 sample - a goal pursued by specific runs of obser- 
vations with Chandra (see V06 and this paper) - we will keep the 
original [O III] line luminosities reported by Z03 and discuss how 
the differences with respect to R08 impact on our conclusions on 
the number density of Compton-thick quasars in §4. 

The outline of the paper is as follows: in §2 the SDSS Type 2 
quasar sample is presented, along with Chandra follow-up observa- 
tions and X-ray spectral results, while §3 is focused on Spitzer ob- 
servations and usage of the mid-IR luminosity as a proxy for the in- 
trinsic X-ray luminosity, hence to pick up heavily obscured quasars. 
The space density of SDSS Compton-thick quasars is shown and 
extensively discussed in §4 in the light of current XRB model ex- 
pectations and recent findings at higher redshifts. Finally, the main 
results are summarized in §5. 

Hereafter we adopt the "concordance" (WMAP) cosmology 
(//o=70 km s-i Mpc"' , aM=0.3, and Ha =0.7; Spergel et al. 2003). 



2 SAMPLE SELECTION AND CHANDRA 
OBSERVATIONS 

The 12 sources of the main sample presented in this work were 
specifically selected for Chandra Cycle 8 observations on the ba- 
sis of their luminous [O III] emission. Combined with previous 
sources targeted by Chandra (see V06 and Ptak et al. 2006), these 
observations provide a complete sampling of the luminous (i.e., in 
the quasar range) Type 2 AGN population revealed by the orig- 
inal study of Z03. In particular, using the M94 correlation be- 
tween [O III] and hard X-ray flux, all the Z03 objects with pre- 
dicted 2-10 keV luminosity above 4 x lO"*"* erg s^', i.e., in the 
quasar regime, have been observed by Chandra and XMM-A'ewton 
(grey region in Fig. [TJ. The "membership" of all these sources 
to the "quasar locus" (above 10^"* erg s^') is valid even adopt- 
ing the lowest X-ray luminosities predicted from the M94 corre- 
lation assuming its lo uncertainty, i.e., the sources included in 
the grey region of Fig. [T] have predicted hard X-ray luminosi- 
ties at least of 10'*'* erg s^'. Only two sources in this "quasar 
locus" were not targeted, because they are radio loud (similarly 
to SDSS 081253.09+401859.9 of the archival sub-sample pre- 
sented here), hence possibly not representative of the majority of 
the Type 2 quasar population, which is expected to be radio quiet. 
The main sample of 12 targets shown in Fig. [T] has been distin- 
guished from the additional three sources retrieved from the Chan- 
dra archive on the basis of the size of the symbols (large and small 
open squares, respectively). 



2.1 Chandra data reduction and analysis 

The 12 SDSS Type 2 quasar candidates of the main sample were 
targeted by Chandra during Cycle 8 with the Advanced CCD Imag- 
ing Spectrometer (ACIS; Garmire et al. 2003) and the S3 CCD at 
the aimpoint; we used a snapshot strategy (~ 10 ks exposures) 
and very faint mode for the event telemetry format. The three 
archival observations, with similar exposures, were carried out dur- 
ing Cycle 7 in faint mode. Standard data reduction procedures 
were adopted using the Chandra Interactive Analysis of Obser- 
vations (CIAO) Version 3.4 software. Source detection in the full 
(0.5-8 keV), soft (0.5-2 keV) and hard (2-8 keV) energy bands 
was carried out with WAVDETECT (Freeman et al. 2002) using a 
false-positive probability threshold of lO^'*, which is less conser- 
vative than that typically adopted in source detection in X-ray sur- 
vey fields, given the "a priori" knowledge of the source positions 
in our case. Limited intervals of high background were present in 
two observations and removed through the generation of new good- 
time intervals. The observation log is shown in Table [T] hereafter 
we will refer to the SDSS sources mostly using their abbreviated 
names. 

In the total sample of 15 sources, only two targets were clearly 
not detected (SDSS J0149-0048 and SDSS J0921+5153), while 
all of the remaining sources were revealed by Chandra, although 
nine with less than 10 counts in the observed full band (see Tz- 
ble[2l(. These objects were carefully checked, following a proce- 
dure similar to that described in Vignali et al. (2001). Monte-Carlo 
simulations indicate that all of the sources with less than 10 counts, 
including the X-ray weak objects with 2~A counts in the detection 
bands, are reliable detections (see last column of Table[2}. From all 
these observations, coupled with those published by V06, we aim 
at deriving the basic X-ray properties (e.g.. X-ray flux, hence lu- 
minosity, for all sources, and possibly photon index and absorption 
for the X-ray brightest targets) of SDSS Type 2 quasars. 
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THE CHANDRA QUEST FOR TYPE 2 QUASARS 

Table 1. Chandra observation log for the SDSS Type 2 quasar candidates presented in this paper 



3 



Src. ID # Object Name z X-ray X-ray Aopi-x Obs. Date Exp. Time 

SDSS J (012000) (82000) (arcsec) (ks) 



Sources observed as targets in Chandra AOS - Main Sample 



9 


005621.72-1-003235.8 


0.484 


00 56 21.7 


+00 32 35.9 


0.2 


2008 Feb 08-09 


9.91 


16 


012032.21-005502.0 


0.601 


01 20 32.2 


-00 55 02.0 


0.2 


2007 Feb 18 


9.83" 


20 


013416.34+001413.6 


0.555 


01 34 16.3 


+00 14 13.7 


0.1 


2007 Sep 10 


9.91 


29 


014932.53-004803.7 


0.566 








2007 Aug 30 


10.01 


30 


015716.92-005304.8 


0.422 


01 57 16.9 


-00 53 04.5 


0.3 


2007 Jun 18 


9.92 


100 


073745.88+402146.5 


0.613 


07 37 45.9 


+40 21 45.7 


0.9 


2007 Feb 03-04 


9.33 


153 


092152.45+515348.1 


0.587 








2007 Sep 27 


10.04 


182 


102746.03+003205.0 


0.614 


10 27 46.0 


+00 32 04.7 


0.4 


2007 Jan 13 


9.91 


186 


103951.49+643004.2 


0.402 


10 39 51.5 


+64 30 04.3 


0.1 


2007 Feb 04 


9.52" 


205 


122845.74+005018.7 


0.575 


12 28 45.7 


+00 50 18.9 


0.2 


2007 Mar 12 


9.42 


232 


144642.29+011303.0 


0.725 


14 46 43.0 


+01 13 03.2 


0.2 


2007 Mar 22 


10.04 


244 


151711.47+033100.2 


0.613 


15 17 11.5 


+03 31 00.3 


0.1 


2007 Mar 28 


9.91 



Archival Sample 



18 012341.47+004435.9 0.399 01 23 41.5 +00 44 35.9 0.2 2006 Feb 07-08 9.83 
117 081253.09+401859.9 0.551 08 12 53.1 +40 19 00.5 0.6 2005 Dec 11 9.89 

152 092014.11+453157.3 0.402 09 20 14.1 +45 31 56.9 0.4 2006 Mar 05 10.05 

The source ID is taken from Table 1 of Z03. The optical positions of the quasars can be drawn from their SDSS name, while 
the X-ray positions for the X-ray detected sources have been obtained from either WAVDETECT or the centroid of the source 
count distribution using the full-band image. We note that also the sources from the archival sample were targets of pointed 
Chandra observations (Cycle 7). " Corrected for high-background periods. 



2.2 X-ray spectral analysis 

Among the 15 targets analyzed in this paper, low-to-moderate 
quality X-ray spectral analysis was possible only for three 
sources, having a number of full-band counts in the range 
fa 50-200: SDSS J1228+0050 (main sample), SDSS J0123+0044 
and SDSS J0812+4018 (archival sample; see TablelU. 

Spectral analysis was carried out with XSPEC Version 11.3.2 
(Amaud 1996) using unbinned data and the C-statistic (Cash 1979; 
Nousek and Shue 1989) for source SDSS J1228+0050, while for 
the other two sources data were binned to at least 10 counts per bin. 
Errors are quoted at the 90 per cent confidence level for one inter- 
esting parameter (AC = 2.71; Avni 1976; Cash 1979), unless stated 
otherwise. Galactic absorption (from Dickey and Lockman 1990; 
see Table[3} was included in all the spectral fittings. 

The limited energy band allowed by Chandra, coupled to 
the poor counting statistics, prevents us from an accurate esti- 
mate of the source spectral parameters. Having said that, we note 
that a fitting with a power-law model provides a flat photon in- 
dex for all of the sources (T = -O.3I+043, T = 
T = 1.22 ±0.22 for SDSS J 1228+0050, SDSS J0123+0044 and 
SDSS J0812+4018, respectively). In the light of the optical clas- 
sification of these sources, the interpretation of the flatness of the 
X-ray spectrum (if compared to that typically expected from Type 1 
quasars; see, e.g., Piconcelli et al. 2005 and references therein) as 
due to absorption at the source provides a reliable representation of 
the X-ray data. Indeed, the inclusion of absorption provides a more 
reliable fitting in term of resulting photon indices, which become 
consistent with those typically observed in Type 1 AGN (Page et 
al. 2OO5)0The derived column densities are 1.52+q^3 x 10^^ cm"^ 
(SDSS J1228+0050), 1.44+i;;l^ x 10^3 cm-2 (SDSS J0123+0044) 



2 The resulting photon index is F = 1.9+[-^ for SDSS J1228+0050 and 



and 2.14+(')-^^ x 10^2 cm^^ (SDSS J0812+4018; see Table[3ll. We 
note that the source with the lowest measured X-ray absorption is 
SDSS J08 12+40 18, which is also radio loud, with a radio loud- 
ness parameteJ3 of ~ 4300. The interpretation of the X-ray data for 
this source is complicated by its radio-loud nature; we note that in 
this case absorption could be ascribed to gas entrained by the un- 
resolved radio jet (which may be responsible also for part of the 
X-ray emission), although no firm conclusions can be drawn from 
the current X-ray data. 

For what concerns the remaining X-ray detected sources, 
we were able to provide some constraints on the absorption 
only for SDSS J0157-0053; its hardness ratio [defined as 
HR=(H— S)/(H+S), where H and S are the source counts in the 
hard and soft bands, respectively] of 0.48±0.32 is suggestive of a 
column density of fa 5 x 10^^ — lO^-* cm^^. 

2.3 X-ray results: discussion 

In this paper, as well as in V06, we have shown that the population 
of SDSS Type 2 quasar candidates investigated by Chandra is char- 
acterized by low X-ray emission, in comparison with that expected 
from the [O III]5007A emission-line intensity and the M94 correla- 
tion. In the following discussion, we consider all of the sources with 
either a Chandra or XMM-Newton observation. We note that for six 
sources for which we refer to Chandra data, XMM-Newton obser- 
vations are also present in the archive; in particular, four of these 
six objects were observed as targets by XMM-Newton. Of these 
six sources, we have only one detection; for the other sources, the 

r = 2.6+15^ for SDSS J0812+4018, while was fixed to F = 2.0 for source 
SDSS J0123+0044. 

' Radio-loudness parameter, defined as i? = /s GHZ//4400 A (where 5 GHz 
and 4400 A are both rest frame; Kellermann et al. 1989). 
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Table 2. X-ray counts, count rates, and results from Monte-Carlo simulations. 



Source X-ray net counts Count rate Det. Sign." 

Name 



SDSSJ [0.5-2 keV] [2-8 keV] [0.5-8 keV] [0.5-8 keV] 





9 n+2.9 
2^-1.6 


< J.U 


9 q+2'9 
^■9-1.6 


z.yj [i.ji— j.ojj X iu 


-4 


3.6 . 


..2.9 


0120—00^5 


1 9+f 5 


1 8+? o 

^■°-I.2 


3 7+fi 

•'^ '-1.8 


3 77 n 93-6 921 x 10" 


-4 


4.0 3.3 4.8 


0134+0014 




<4.8 


2 7+2'' 

^■'-1.5 


2.73 [1.21-5.85] xlO" 


-4 


3.6 . 


..4.0 


0149-0048 


<4.8 


< 3.0 


<4.8 


<4.80 xlO-'* 








0157-0053 






99 7+5.8 


2.29 [1.82-2.87] xlO" 


-3 






0737+4021 


1 Q+26 
^■^-1.3 




4 7+3.3 
^■'-2.1 


5.03 [3.00-8.58] x 10" 


-4 


3.7 4. 


1 5.5 


0921+5153 


<4.8 


< 3.0 


<4.8 


<4.78 xlO-** 








1027+0032 


2 9+2' 


< 3.0 


2 7+2'' 

^■'-1.5 


2.72 [1.21-5.85] xlO" 


-4 


4.7 . 


. . 3.3 


1039+6430 




<6.4 


•S n+3.6 


6.20 [3.68-9.98] x 10" 


-4 


5.7 . 


..5.7 


1228+0050 






50.4+^-^ 


5.35 [4.60-6.22] x 10" 


-3 






1446+0113 


5 9+3 'i 


< 3.0 


5.8li;f 


5.78 [3.49-9.37] x 10" 


-4 


6.4 . 


..5.5 


1517+0331 


2 9+-' 


< 3.0 


^■vli;? 


2.72 [1.21-5.65] xlO" 


-4 


4.4 . 


. . 3.8 


0123+0044 


13 5+* '* 

"■■'-3.6 


145.2^!3;> 


i6i.2;i3-? 


1.64 [1.51-1.78] xlO- 


-2 






0812+4018 


128.6+1?'^ 


72.0^^ 


199.8115-^ 


2.02 [1.88-2.17] xlO- 


-2 






0920+4531 




< 6.4 


Q 4+4.2 
^■^-3.0 


9.36 [6.37-13.54] xlO 


-4 


7.2 . 


..7.9 



The source names are abbreviated on the basis of the four RA and DEC digits. EiTors on the X-ray net 
counts (i.e., background-subtracted) were computed according to Gehrels (1986). The upper limits are at 
the 95% confidence level and were computed according to Kraft, Burrows & Nousek (1991). " Significance 
of the detection (in units of o) in the soft band, hard band, and full band, respectively, estimated using the 
Monte-Carlo method described in Vignali et al. (2001) and applied only to sources with less than 10 counts. 



Table 3. Propeilies of the SDSS Type 2 quasars observed by Chandra. 



Object 




log ^0l^] 


^2-10 keV 




^2-10 keV 


Ll-utkeV (pr.) 


Si A GHz 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) (9) 


0056+0032 


2.79 


9.45 


1.0x10-'^ 




8.9x10"' 


1.5xlO""-2.6xlO"5 


8.60 8.33 


0120-0055 


3.69 


9.28 


1.4x10-15 




2.1xlO« 


1.0xlO""-1.8xlO"5 


< 0.144 


0134+0014 


2.48 


9.53 


9.6x10-'*' 




1.2xlO''2 


1.8xlO""-3.2xlO"5 


< 0.201 


0149-0048 


2.55 


9.36 


< 1.7 X 10-'^ 




< 2.2 X 10"2 


1.2xlO""-2.1xlO"5 


1.03 1.60 


0157-0053 


2.58 


9.52 


8.1x10-'^ 




5.2xlO"2 


1.8xlO""-3.1xlO"' 


< 0.135 


0737+4021 


6.18 


9.28 


1.9x10-'"^ 




3.0xlO*'2 


1.0xlO""-1.8xlO"5 


< 0.134 


0921+5153 


1.27 


9.28 


< 1.6 X 10-'^ 




< 2.3 X 10"2 


1.0xlO""-1.8xlO"5 


2.37 2.49 


1027+0032 


4.47 


9.36 


1.0x10-'^ 




1.6xlO"2 


1.2xl0""-2.1xl0"5 


3.25 3.09 


1039+6430 


1.18 


9.41 


2.1x10-'^ 




1.2xl0"2 


1.4xl0""-2.4xl0"5 


<0.5 


1228+0050 


1.88 


9.28 


1.8x10-'* 


1.52^i;30 X 10^3 


3.5xlO''3 


1.0xlO""-1.8xlO"5 


3.14 2.95 


1446+0113 


3.55 


9.37 


2.1x10-'^ 


5.0xlO"2 


1.3xlO""-2.2xlO"5 


5.75 5.75 


1517+0331 


3.78 


9.36 


9.8x10-'* 




1.5xlO"2 


1.2xl0""-2.1xl0"5 


< 0.150 


0123+0044 


3.24 


9.13 


5.8x10-'" 


IMtlf, X 10^3 


3.4x10"" 


7.3xlO"3-1.3xlO"5 


11.83 11.34 


0812+4018 


5.16 


9.39 


7.5x10-'" 


2A4toil X 10^2 


1.7x10"" 


1.3xlO""-2.3xlO"5 


1067.88 1058.83 


0920+4531 


1.51 


9.04 


3.2x10-'^ 


1.8xlO"2 


5.9xlO"3-1.0xlO"5 


<0.5 



(1) Abbreviated SDSS name; (2) Galactic column density, from Dickey & Lockman (1990), in units of 10^" cm-^; (3) log of the [Olll] line 
luminosity, in units of L ;, (from Z03); (4) Galactic absorption-corrected flux (or upper limit, in units of erg cm-" s-') in the 2-10 keV band, 
extrapolated from the observed 0.5-8 keV count rate or upper limit assuming a power law with T = 2.0 (typical for AGN X-ray emission). We 
note that assuming a pure cold reflection spectrum, parameterized here as a power-law with F = 0, would increase the observed 2—10 keV flux 
by a factor of ft; 8; (5) intrinsic column density, obtained from the best-fit X-ray spectral fitting, when possible, in units of cm-^; (6) 2-10 keV 
rest-frame luminosity, obtained through the observed flux (or upper limit), connected for the effect of intrinsic absorption (when the column 
density can be measured directly from the spectral fit), in units of erg s-' ; in case of a reflection spectrum, the rest-frame 2-10 keV luminosity 
would increase by a factor of Ri 3. (7) 2-10 keV luminosity range predicted from the [O III] line vs. hard X-ray flux correlation (M94), in units 
of erg S-' ; the interval is given by the ±la values around the mean M94 correlation; (8)-(9) integrated (8) and peak (9) radio flux density (or 
l0 upper hmits) at 1.4 GHz, all from the FIRST (Becker, White & Helfand 1995) except for SDSS J1039+6430, whose upper limit comes 
from the NVSS (Condon et al. 1998), in units of mJy. We note that for SDSS J1027+0032 the optical-radio association is dubious, since the 
distance between the radio and optical position is large (1.0"). 
All luminosities are computed using Ho=70 km S-' Mpc-' , nM=0.3, and 0.^=0.1. 



© 2010 RAS, MNRAS 000,[T]-?? 



THE CHANDRA QUEST FOR TYPE 2 QUASARS 5 



4x 



A 

H 



® 



• * • >' ML* 



□ This work 
/\ V06 sample 
O V04 sain pic 



- lO*' 



- 10« 



0.4 0.6 
Rcdshift 



0.8 



1046 



Si 045 



Corripton thick? 




H This work (spec, analysis)- 
Q This work (no spectra) 
Av06 (spec, analysis) - 
Z^V06 (no spectra) 



10« 



10« 



10« 



1045 



1046 



^-10 kcv 



(mcasurcd/dc absorbed) 



Figure 1. Logarithm of the measured Lp in| luminosity (not corrected for 
extinction within the narrow-line region) vs. redshift for all of the sources 
in the Zakamska et al. (2003) catalog (small filled circles). At the right of 
the panel, the 2-10 keV luminosity, estimated using the [O III] emission 
and the M94 correlation, are shown. The key provides a description of the 
X-ray observations; in particular, the sources presented in this paper are 
plotted as open squares (the three small open squares refer to the objects 
whose observations were retrieved from the Chandra archive). The grey 
region defines the locus of the AGN with log Lx> 44.6 erg s"', where the 
AGN still he in the quasar luminosity range (log Lx> 44 erg s^') even 
assuming the la error around the mean of the original correlation between 
[O III] and hard X-rays. We note that all of the objects in the grey region 
have been observed in X-rays except for two sources which are radio loud. 

integration time is typically lower or comparable to that of Chan- 
dra observations, once the intervals of high particle background 
(flares) are removed. This means that the constraints achieved us- 
ing Chandra data for these six sources remain the best available to 
date. About the source detected by XMM-Newton, it was already 
detected by Chandra (src. #15; V06) and was bright enough for 
an X-ray spectral extraction. However, we note that the presence of 
diffuse X-ray emission in the XMM-A'ewton observation within the 
aperture adopted for X-ray spectral extraction may provide some 
contamination to the AGN properties, as shown by the source flux 
which is ~ 70 per cent higher than in the Chandra observation. For 
this reason, we prefer to use the Chandra constraints (as reported 
in V06) also for this source. 

If we correct the X-ray luminosity accounting for the mea- 
sured absorption (when a column density estimate is feasible 
through spectral analysis, see the filled squares and triangles in 
Fig. |2j the measured column densities are in the range ~ 1.5 x 
10^2 - 4.3 X 10^3 cm-2, see V06 and ^n\ . we generally find a 
good agreement between the expected and the observed 2-10 keV 
luminosity. On the other hand, there is a significant fraction of 
Type 2 AGN candidates (~ 65 per cent) whose observed luminos- 
ity is less than 1 per cent of the predicted one; for these sources, 
heavy obscuration towards the nuclear source appears a reason- 
able explanation for their weak X-ray emission, supported by the 
optical classification of Z03 (see also Ptak et al. 2006 for similar 
conclusions, and the recent results of LaMassa et al. 2009 using a 



Figure 2. Comparison of the 2-10 keV rest-frame luminosity coinputed 
from the available X-ray data with that predicted assuming the M94 cor- 
relation (i.e., derived from the [O III] luminosity). The dashed line shows 
the 1:1 ratio between the two X-ray luminosities. All of the sources with 
either a Chandra or XMM-A'ewton observation are shown. For the filled 
symbols (squares: data presented in this paper; triangles: data in V06), the 
X-ray luminosity is de-absorbed on the basis of the column density de- 
lived from the best-fitting model (see Table [5); for all of the remaining 
sources, the X-ray luminosity is derived from the X-ray flux with no cor- 
rection for the unknown column density. Sources with X-ray spectral results 
at the center of the plot are obscured, with typical measured column densi- 
ties of a few X 10^2-5 x 10^3 cm ^, while the upper-left grey region shows 
the locus of extremely obscured, likely Compton-thick sources, where the 
observed luminosity is <1 per cent of the predicted one. Left-ward arrows 
indicate upper limits on the observed X-ray luminosity. 

lower redshift sample of SDSS Seyfert 2 galaxies). Hereafter, we 
will assume the sources with measured/predicted X-ray luminosity 
less than 0.01 as probably Compton thick. This is a conservative 
threshold since there are known Compton-thick AGN in the local 
Universe with more than 1 per cent of observed over intrinsic X-ray 
emission (e.g., NGC 6240 and the Circinus galaxy; Bassani et al. 
1999; Matt et al. 2000; see also Maiolino et al. 1998, where this 
fraction is ~ 2 per cent for a sample of Compton-thick Seyfert 2 
galaxies). 

Similarly to V06, we did not correct the [O III] flux to ac- 
count for the absorption due the narrow-line region (NLR) itself 
(for details, see Maiolino et al. 1998; Bassani et al. 1999; Panessa 
et al. 2006, hereafter P06; see also Lamastra et al. 2009), because 
the SDSS spectral coverage does not typically cover the rest-frame 
wavelengths necessary to measure the emission-line fluxes of HP 
and Ha, hence to use the Ha/HP ratio as extinction measurement. 
Using the ratio of higher order Balmer lines (e.g., Hp/Hy, Hy/HB) to 
estimate extinction represents a difficult task, due to their weakness 
over often noisy spectra. This approach implies that the predicted 
X-ray luminosities would be even higher for at least some sources, 
thus generally requiring a larger column density to account for the 
observed luminosities. 

Furthermore, we have assumed the correlation between the 
[O III] emission-line intensity and the hard X-ray emission from 
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M94, which is based on a large sample (> 100) of local Seyfert 
galaxies. As reported in TablejS] our objects have X-ray luminosi- 
ties well below those predicted by the M94 correlation; 13 out 
of the 15 targets have X-ray measured/predicted luminosity ratios 
(Lx^meas/Lx [oill]) '^e range ~ 0.001 — 0.081 (where for the pre- 
dicted X-ray luminosity, the "central" value from the M94 correla- 
tion has been assumed). 

A few caveats must be considered before comparing different 
[O III] vs. hard X-ray correlations. If a pure cold reflection spec- 
trum, roughly parameterized here as a flat (F = 0) power-law slope, 
were present in the observed 0.5-8 keV band covered by Chan- 
dra, the 2-10 keV fluxes reported in Table[3]would increase by a 
factor ~ 8, and the rest-frame 2-10 keV luminosities by a factor 
3. Still a significant (fa 40-50 per cent) fraction of SDSS Type 2 
quasars would populate the Compton-thick locus in Fig.[2] We note, 
however, that studies of local reflection-dominated AGN indicate 
that a scattering component is possibly present in such cases (e.g., 
Comastri 2004), mostly "contaminating" the soft X-ray band; this 
would produce softer spectra and, as a consequence, any column 
density estimated through low-to-moderate quality X-ray spectral 
information would be under-estimated. 

Since M94 work, many authors have searched for the presence 
of a correlation between the [O III] emission-line intensity and the 
hard X-rays. H05 clearly show the differences in deriving such cor- 
relation from an optically selected sample. In particular, they found 
that log(L2-i0i(:ev/i[o III]) 1-59 (dispersion 0=0.48) for optically 
([Olll]) selected Type 1 AGN, while log(L2_ioi:fv/^[o iii])=0-57 
(dispersion 0=1.06) for optically selected Type 2 AGN. We re- 
mind that the best-fit ratio assumed throughout this paper is 
^2-l0jt(>v/i-[O lll]~ 1-8 (o ~ 0.6), similar for Type 1 and Type 2 
AGN (M94); this value is not significantly different from that of 
Type 1 AGN in H05. Under the basic assumption (unified schemes) 
that the intrinsic i2-l0*:cv/^[o ill] should be the same for Type 1 and 
Type 2 AGN, the only difference being the obscuration/extinction 
of the nuclear emission, to report our predicted X-ray luminosities 
(reported in Table [3] and shown in Fig. |2} into those of H05 we 
have to divide them by ~ 1.5. Given the significant dispersion of 
H05 correlation, a large fraction of sources (~ 50 per cent) would 
still be in the locus of heavily obscured, possibly Compton-thick 
quasars. 

Finally, we note that the correlation found by P06 produces 
even larger hard X-ray luminosities than those reported in this pa- 
per (see also For consistency with the assumptions of V04 and 
V06 which, combined with the present work, provide a complete 
X-ray coverage of all radio-quiet Type 2 quasars with logLjoni] > 
9.28 Lq selected from Z03, we have decided to adopt the M94 cor- 
relation, keeping in minds all the caveats concerning both the selec- 
tion criteria of the different samples and the significant dispersion 
in all of the correlations presented above; for further discussion, re- 
fer to ^ Support to the reliability of the [O IIl]-X-ray correlation 
is also provided by Spitzer data (see 



3 NEAR AND MID-IR OBSERVATIONS OF SDSS TYPE 2 
QUASARS 

3.1 Spitzer data reduction and analysis 

In the following, we estimate the nuclear power of the AGN by 
means of the rest-frame 5.8nm and 12.3nm luminosities, providing 
a robust support to the results from X-ray analysis. To this purpose, 
we retrieved from the Spitzer archive both IRAC and MIPS photo- 
metric observations for 20 SDSS Type 2 quasars in common with 



those reported in this work and in V04a and V06. Most of these 
sources were targeted by Spitzer to get insights, for the first time, 
on the near-IR and mid-IR properties of a well defined sample of 
optically selected Type 2 AGN. 

In particular, by these infrared data, we intend to provide an 
adequate complement to the X-ray investigation presented in the 
previous sections. The X-ray properties of these sources have been 
published mostly in V06 and Ptak et al. (2006); for seven of these 
sources. X-ray data are presented here. While IRAC covers the 
near-infrared (near-IR) and mid-IR bands with four channels (3.6, 
4.5, 5.8 and 8|i.m), MIPS allows for investigations at longer wave- 
lengths (24, 70 and 160nm), although in this work only 24(i.m data 
have been used. 

For both IRAC and MIPS, we used the final post-calibrated 
data (PBCD) produced by the Spitzer Science Center (SSC) 
pipeline (Version 14.0.0 and 16.1.0, respectively). The flux den- 
sities in the IRAC and MIPS bands were measured on the signal 
maps using aperture photometry. The chosen aperture radius for 
the IRAC bands was 2.45", and the adopted factors for the aperture 
corrections are 1.21, 1.23, 1.38 and 1.58 (following the IRAC Data 
Handbook) at 3.6nm, 4.5iim, 5.8nm, and 8nm, respectively. The 
aperture radius for the 24|i.m MIPS band was 7", and the resulting 
aperture flux density was corrected adopting a factor of 1.61 (see 
the MIPS Data Handbook). 

All of the quasars have been detected by Spitzer, spanning a 
flux density range at Syum (24(j.m) of ~ 0.11 mjy (~ 1.46 mjy; 
source ID=290, the faintest object in the current sample) to 
~ 12.8 mJy (~ 56.4 mJy; source ID=197, the brightest source). 
Only four sources show indications of limited contamination by 
a nearby companion; however, in all of these cases the nearby 
source shows a decreasing flux density at longer wavelengths, sug- 
gesting a stellar or galaxy origin for its near-IR/mid-IR emission. 
The rest-frame 5.8nm and 12.3ti.m luminosities (Table|4]l were de- 
rived via simple interpolation between nearby bands in the source 
rest frame; this method avoids that the results are dependent upon 
the assumption of particular AGN (torus) templates to reproduce 
the observed data points. Furthermore, we checked our results 
against the 14.5(i.m luminosities obtained using IRS spectra and 
reported by Zakamska et al. (2008) for eight sources in common, 
and found a relatively good agreement. The combined statistical 
and systematic errors associated to these mid-IR luminosities are 
~ 10-15 per cent. In Table]?] also basic X-ray information and cor- 
responding reference paper are presented. 



3.2 Mid-IR luminosity as a proxy of the X-ray emission 

The rest-frame 5.8iim and 12.3|i.m luminosities (reported in Ta- 
ble|4]l have been used to estimate the intrinsic, rest-frame 2-10 keV 
luminosities (Lx,mid-IR) assuming the relations provided by F09 
and Gandhi et al. (2009; hereafter G09), respectively. 

The former relation^ parameterized by logL2_i()keV = 
43.57+0.72x(logl58-44.2), has been calibrated using Type 1 
AGN in the Chandra Deep Field-South (CDF-S; Brusa et al. 
2010) and in the Chandra Cosmic Evolution Survey (C-COSMOS; 
Civano et al., in preparation); for further details, see §3.3 and Fig. 5 
in F09. This relation is similar (both in slope and intercept) to 
that recently derived by Lanzuisi et al. (2009) from a sample of 



* This relation is valid if the criterion log 158>43.04 is satisfied, as in the 
case of all of our sources; see F09 for details. 



© 2010 RAS, MNRAS 000,[T]-?? 



THE CHANDRA QUEST FOR TYPE 2 QUASARS 7 



Table 4. Properties of the SDSS Type 2 quasars observed by Spitzer IRAC and MIPS. 


ID 


Src. Name 


z 


log 158 


Ci V J I L.J 


X-ray into 


log Lx 








iLy,) 


(Lq) 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


9 


005621.72+003235.8 


0.484 


11.3 


11.6 


D/CA'09 


41.95 


18 


012341.47+004435.9 


0.399 


11.2 


11.2 


D+/CA'09 


44.53 


20 


013416.34+001413.6 


0.555 


11.2 


11.5 


D/C/V09 


42.08 


30 


015716.92-005304.8 


0.422 


10.7 


10.9 


D/CA'09 


42.72 


34 


021047.01-100152.9 


0.540 


11.2 


11.3 


D/XA'Oe 


44.34 


83 


031950.54-005850.6 


0.626 


11.0 


11.2 


D/CA'Oe 


42.57 


113 


080154.24+441234.0 


0.556 


11.2 


11.4 


D+/C/V06 


44.62 


117 


081253.09+401859.9 


0.551 


11.4 


11.4 


D+/CA'09 


44.23 


119 


081507.42+430427.2 


0.510 


11.8 


11.8 


ND/CA'06 


< 42.51 


130 


084234.94+362503.1 


0.561 


10.9 


11.4 


ND/CA'06 


< 42.89 




00901 4 1 1 -l-ilS^ 1 S7 ^ 


402 


1 1 3 


1 1 s 

1 L.J 




42 26 


186 


103951.49+643004.2 


0.402 


11.6 


11.7 


D/CA'09 


42.08 


196 


115314.36+032658.6 


0.575 


11.2 


11.3 


D-I-/C/V06 


44.08 


197 


115718.35+600345.6 


0.491 


12.1 


12.2 


ND/C/V06 


< 42.18 


204 


122656.48+013124.3 


0.732 


11.6 


11.6 


D-I-/X/V06 


44.63 


207 


123215.81+020610.0 


0.480 


11.1 


11.3 


D/C/V06 


42.30 


239 


150117.96+545518.3 


0.338 


10.7 


10.9 


D/R/V04 


44.24 


256 


164131.73+385840.9 


0.596 


11.6 


11.8 


D+/X/WQ6 


44.87 


289 


235818.87-000919.5 


0.402 


10.6 


10.9 


ND/C/V06 


< 42.49 


290 


235831.16-002226.5 


0.628 


10.3 


10.8 


D/C/V06 


42.51 



(1) Source ID (from Table 1 of Z03); (2) SDSS source name; (3) source redshift; (4) log of the luminosity at 
5.8 i^m (rest frame), in units of Lv,; (5) log of the luminosity at 12.3 i^m (rest frame), in units of L :,; (6) avail- 
able X-ray information: the first term indicates whether the source has been detected in X-rays (D: detected; 
D+: detected and X-ray spectrum available; ND: undetected); the second term indicates the X-ray satellite 
whose data have been used to obtain information (C: Chandra; X: XMM-Newton; R: i?05Ar/RASS). The 
last term indicates the reference paper (where V09 means the current work); (7) logarithm of the rest-frame 
2-10 keV luminosity (erg s^') or upper limit when the source is undetected. When an X-ray spectrum is 
available, the X-ray luminosity has been corrected to account for the measured absorption. 
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Figure 3. Comparison of the 2-10 keV luminosity computed from the available X-ray data with that predicted from the mid-IR luminosity at 5.8i^m (rest 
frame) assuming the Fiore et al. (2009) congelation {left panel) and at 12.3iJ.m (rest frame), using the Gandhi et al. (2009) correlation {right panel). The dotted, 
short-dashed and long-dashed lines indicate ratios of 1:1, 1:10 and 1:100 between the measured and the predicted X-ray luminosity. Symbols are the same as in 
Fig. [2] the filled circle marks the source with ROSAT (RASS) detection (V04). The letter inside each symbol indicates whether the source is a Compton-thick 
(T) or a Compton-thin (t) AGN candidate on the basis of the [O III]-driven results presented in j|2.3l and shown in Fig. 121 for one object, the analysis does not 
lean towai'ds any of the two classifications [hence we have an "intermediate" (i) AGN], although HR analysis suggests Compton-thin absorption. 
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SWIRE sources with X-ray coverage. For comparison, the assump- 
tion of the Lutz et al. (2004) relation between the 6iim and hard 
X-ray luminosity derived using a sample of 71 local, moderate- 
luminosity Seyfert galaxies with low-resolution ISO spectra would 
predict larger X-ray luminosities given the mid-IR luminosities of 
our sample (see below for details). 

The latter relation, parameterized by logL2_iokeV = 
(log(A.Lx)i2.3nm+4.37)/l.l 1, has been obtained from a sample of 42 
Seyfert galaxies with near diffraction-limited mid-IR imaging ob- 
tained with the VISIR instrument at the VLr(G09; see also Horst 
et al. 2006, 2008 and 2009). 

All of these correlations are based on the assumption that the 
mid-IR band provides a clear view of the AGN nuclear activity, 
which is more affected by extinction and absorption at optical, ul- 
traviolet and X-ray frequencies. However, a few issues must be 
taken into account in the following discussion. X-ray information 
is typically obtained from all the available data in literature work 
or in archives (from ASCA to BeppoSAX, up to the current X-ray 
missions); we note, however, that F09 and Lanzuisi et al. (2009) 
use only data with moderate (a few arcsec; XMM-A'evvton) to high 
(below 1"; Chandra) resolution, providing overall a relatively good 
match with Spitzer imaging data in terms of angular resolution|£|On 
the one hand, G09 mid-IR data are characterized by sub-arcsec res- 
olution (corresponding to physical scales of less than half parsec) 
matched to a broad variety of X-ray data, including old {ASCA) 
and recent (Suzaku) observations. While host galaxy contribution 
to hard X-ray emission, at the AGN luminosities probed by most 
of their sources (above lO'*^ erg s^'), is likely negligible (although 
some reprocessed emission may be present), it can still contaminate 
the mid-IR data at some level. On the other hand, at the median red- 
shift (z fa 1.5) of the Type 1 AGN used by F09, the mid-IR emis- 
sion would probe physical scales of ~ 20 kpc. Finally, Lutz et al. 
(2004) ISO spectra are characterized by 24" x 24" apertures (corre- 
sponding to physical scales of a few kpc), including, in many cases, 
noticeable host galaxy emission, as pointed out by the authors. 

The results are shown in Fig. [3] where the 2-10 keV lumi- 
nosity, predicted using the luminosity at 5.8|J.m (left panel) and 
12. 3 Jim (right panel) is plotted against the measured X-ray lumi- 
nosity; with the term "measured X-ray luminosity" (similarly to 
Fig. [2}, we mean that the source luminosity has been corrected for 
the effect of absorption, once X-ray spectral analysis is able to pro- 
vide a reliable estimate for it (filled squares: this work; filled trian- 
gles: V06); in all of the remaining cases, no correction has been ap- 
plied (open squares and triangles). In addition, the filled circle rep- 
resents source ID 239, which was detected by the ROSAT All Sky 
Survey (RASS; V04). The dotted, short-dashed and long-dashed 
lines indicate ratios of 1:1, 1:10 and 1:100 between the measured 
and the predicted X-ray luminosity. The average X-ray luminosi- 
ties predicted by the two methods are log Lx=44.02 erg s^' and 
44.48 erg s^' using the F09 and G09 mid-IR vs. X-ray correla- 
tion, respectively. For comparison, the average X-ray luminosity 
assuming the Lutz et al. (2004) relation is log Lx=44.41 erg s^', 
15 per cent lower than that predicted using G09 correlation. 

Overall, it appears evident that in both plots two main source 
groups are present: the right-most one comprises sources which are 
more or less consistent with the 1:1 correlation, especially using 
the (A.Lx)i2.3|im as a proxy of the nuclear emission (see below for 
an extended discussion). X-ray spectral analysis indicates that these 



^ We note that source extraction regions of a few tens arcsec are typically 
used in XMM-Newton observations. 



sources are Compton-thin AGN and are the same residing close to 
the 1:1 correlation in Fig. [2] The left-most group regards sources 
where the expected 2-10 keV luminosity is a factor ~ 10-100 
higher than the measured X-ray luminosity. All but one of these 
sources are also among the most extreme objects (i.e., candidate 
Compton-thick AGN) shown in Fig.|2j the only exception is source 
ID 30, which seems to be Compton thin from hardness-ratio analy- 
sis (see 32. 2t . However, it is evident from Fig.[3](left panel) that, for 
the Compton-thin Type 2 quasars of our sample, using the correla- 
tion proposed by F09 provides an offset by a factor ~ 3 with respect 
to the X-ray measured/intrinsic 1 : 1 luminosity ratio, while the G09 
correlation seems to produce a good match with the luminosity ra- 
tio of unity (Fig. [3] right panel), hence with the results obtained 
from [O III] analysis. Although clearly mid-IR vs. X-ray correla- 
tions need to be investigated further on larger samples to achieve a 
more accurate "calibration", we note that, at the measured X-ray lu- 
minosities sampled by the SDSS Compton-thin Type 2 quasars (i.e., 
above lO'*^ erg s^'), the dispersion in the F09 correlation seems to 
be large (see their Fig. 5). By contrast, the correlation presented by 
G09 has a relative small dispersion (0.36 dex) throughout the en- 
tire luminosity range, being smaller for the subsample of 22 well 
resolved sources (0.23 dex; see §4.1 in G09 for details). Therefore, 
in the following section, we will report on the results obtained using 
the information from the [O III] and 12.3vtm to predict 2-10 keV 
luminosities. 

Caution is obviously needed when dealing with infrared data, 
given the likely presence of both AGN (torus) and starburst emis- 
sion. Although the AGN seems to dominate over the Spitzer bands 
(eight of the 20 sources listed in Table |4] have IRS coverage; see 
§3.4 of Zakamska et al. 2008 for details), we note that in the de- 
rived rest-frame 5.8|J.m and 12.3vtm luminosities there might be 
some contribution from star formation (see Table 1 of Zakamska et 
al. 2008). However, we remind that the correlation reported by G09 
was obtained using near diffraction-limited mid-IR imaging, where 
any non-AGN contribution can be observed and kept under control 
(see also Horst et al. 2009). 

Overall, both [O III] and 12.3iim luminosity provide similar 
indications for what concerns the X-ray emission of SDSS Type 2 
quasars, with roughly half of the sampled (i.e., with X-ray cover- 
age) population being possibly Compton-thick (see ^for further 
comparison and discussion). Recently, heavy obscuration has also 
been found in some SDSS Type 2 AGN by Lamastra et al. (2009) 
using XMM-Newton observations; given the low redshifts of their 
sources, they were able to correct for the extinction within the NLR 
and provide a tight and robust [O III] vs. hard X-ray correlation. 

Finally, we note that Spitzer (IRS) follow-up observations 
of SDSS Type 2 quasars have recently revealed a wealth of 
mid-infrared properties, with 7/12 sources showing a silicate ab- 
sorption features at 10ii.m, suggestive of heavy obscuration to- 
ward the source (Zakamska et al. 2008). It is intriguing that the 
two objects with the most prominent silicate absorption features 
(SDSS 0056+0032 and SDSS 0815+4304; see Fig. 1 of Zakamska 
et al. 2008) are candidate Compton-thick quasars accordingly to all 
methods adopted in this paper to estimate the intrinsic X-ray lu- 
minosity, although not all of the SDSS Type 2 quasars which have 
been classified as Compton-thick according to our analyses were 
found to show Si absorption (see Fig. 1 of Zakamska et al. 2008). 
This fact can be explained in clumpy torus models (e.g., Nenkova, 
Ivezic & Elitzur 2002; Honig, Beckert & Ohnaka 2006; Nenkova et 
al. 2008a,b), where the effective infrared optical depth can be de- 
creased by random mis-alignments and holes through the obscuring 
clouds. 
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Table 5. Classification of Type 2 quasars with logLpm] > 9.28 L©. 


ID 


Src. Name 




^X.meas/'^X.fOIII] 


^ A , m e a b ' A . J i . .1 (0. ni 


Class. [O III] 


Class. 12.3|J.ni 






\^) 


('41 




(fi) 


IT] 


1 


005009.81-003900.6 




0.37 




thin 




9 


005621.72+003235.8 


0.484 


1.41E-03 


7.41E-04 


THICK 


THICK 


16 


012032.21-005502.0 


0.601 


4.90E-03 




THICK 




20 


013416.34+001413.6 


0.555 


1.58E-03 


3.39E-03 


THICK 


THICK 


29 


014932.53-004803.7 


0.566 


<4.27E-03 




THICK 




30 


015716.92-005304.8 


0.422 


7.08E-03 


5.01E-02 


interm 


thin 


34 


021047.01-100152.9 


0.540 


0.16 


0.87 


thin 


thin 


83 


031950.54-005850.6 


0.626 


4.37E-03 


1.78E-02 


THICK 


thin 


100 


073745.88+402146.5 


0.613 


7.08E-03 




THICK 




113 


080154.24+441234.0 


0.556 


0.49 


1.3 


thin 


thin 


119 


081507.42+430427.2 


0.510 


<3.89E-03 


<4.79E-03 


THICK 


THICK 


130 


084234.94+362503.1 


0.561 


<2.75E-03 


<2.88E-02 


THICK 


interm 


153 


092152.45+515348.1 


0.587 


<5.37E-03 




THICK 




182 


102746.03+003205.0 


0.614 


3.09E-03 




THICK 




186 


103951.49+643004.2 


0.402 


2.09E-03 


2.00E-03 


THICK 


THICK 


196 


1 15314.36+032658.6 


0.575 


0.13 


0.47 


thin 


thin 


197 


115718.35+600345.6 


0.491 


<1.70E-03 


<1.02E-03 


THICK 


THICK 


204 


122656.48+013124.3 


0.732 


0.42 


0.89 


thin 


thin 


205 


122845.74+005018.7 


0.575 


8.13E-02 




thin 




207 


123215.81+020610.0 


0.480 


1.82E-03 


8.51E-03 


THICK 


THICK 


232 


144642.29+011303.0 


0.725 


9.55E-03 




THICK 




244 


151711.47+033100.2 


0.613 


2.95E-03 




THICK 




256 


164131.73+385840.9 


0.596 


0.40 


1.0 


thin 


thin 


289 


235818.87-000919.5 


0.402 


<6.61E-03 


<2.88E-02 


THICK 


interm 


290 


235831.16-002226.5 


0.628 


1.58E-03 


3.98E-02 


THICK 


thin 



(1) Source ID (from Table 1 of Z03); (2) SDSS source name; (3) source redshift; (4) ratio of the measured 2-10 keV 
luminosity and that predicted from the [O III]5007A luminosity assuming the M94 correlation; (5) ratio of the 
measured 2-10 keV luminosity and that predicted from the 12.3^m luminosity (when available; see Table 
assuming the G09 correlation; (6) tentative source classification on the basis of the measured/predicted (from the 
[O III]) luminosity ratio. For source #30, although a Compton-thick classification is derived from the luminosity 
ratio, hardness-ratio analysis ( j|2.2t suggests Compton-thin obscuration; (7) tentative source classification on the 
basis of the measured/predicted (from the 12.3jim) luminosity ratio; objects are classified as "intermediate" if the 
upper limit on the luminosity ratio prevents any conclusion. 



4 NUMBER DENSITY OF SDSS COMPTON-THICK 
QUASARS 

It is now interesting to compute the space density of candidate 
Compton-thick Type 2 quasars selected by the SDSS and compare 
it with that of other samples of Compton-thick AGN candidates and 
with the predictions by synthesis models of the XRB. First, we con- 
sidered the complete sub-sample of radio-quiet objects presented in 
V06 and in this work with logLjo iu]> 9.28 (therefore source 
SDSS 0812+4018 has been removed because it is radio loud; 
see six objects fall in the redshift range z = 0.3 — 0.5 [with 

5 objects being Compton-thick candidates using the X-ray mea- 
sured/predicted (from [O III]) luminosity ratio, Lx.meas/i'X [Olll] 1 
and 19 objects fall in the redshift range z = 0.5 — 0.8 (with 12 ob- 
jects being Compton-thick candidates, still according to [O III]); 
the properties of these sources are summarized in Table |5] 

Before entering into the details of the number density com- 
putation, we remind the reader that we refer to a Compton-thick 
quasar when the ratio of the measured to predicted X-ray lumi- 
nosity is less than 0.01. On the basis of the ix.meas^^'X [Olll] 
tio, the X-ray emission for 17 of the 25 Type 2 quasars in our 
sample (68 per cent) appears to be obscured by Compton-thick 
material. However, since the Type 2 quasar sample is selected 
on [O III] luminosity, the estimated Compton-thick quasar frac- 
tion may be overestimated. We can provide an indication of the 



"[O III] bias" by determining the fraction of the Type 2 quasars 
that are classified as Compton thick based on both the X-ray/[0 III] 
and X-ray/mid-IR ratios. Of the 16 Type 2 quasars with archival 
Spitzer observations and logLjg ni] > 9-28 Lc-, (see Table [5]l, 10 
have Z-x.meas/^x [Olll] <0.01, of which 6 (60 per cent) also have 
^x,meas/ix,mid-lR<0.01 (where mid-IR is based on 12.3[im); how- 
ever, two objects have X-ray/mid-IR upper limits which are slightly 
above our Compton-thick quasar definition ("intermediate" in Ta- 
ble |5j, suggesting that the fraction of Compton-thick quasars may 
be 75 per cent (i.e., 6 out of 8 objects). 

Overall, our analyses therefore suggest that the unbiased 
Compton-thick quasar fraction is ~ 50 per cent0 

In our calculation below, we will refer to the [O III] indicator 
to compute the number density of SDSS Compton-thick quasars. 



* The current study indicates that, for the 14 sources with [O III] and 
12.3|J.m indicators (16 objects minus two classified as "intermediate" in the 
mid-IR), [O III] is highly reliable in classifying sources as Compton-thick 
in at least 75 per cent of the cases (since source #83 has a "borderline" clas- 
sification as Compton thin according to the 12.3|xm indicator). Applying 
this "coiTection factor" to the number of Compton-thick sources with only 
[O III] as reliable indicator (i.e., no solid mid-IR information) implies at 
least 7 Compton-thick quasars out of 9 (hence, 13/25 for the sample under 
investigation). 
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redshift 

Figure 4. Volume density of Compton-thick AGN samples with different 
intrinsic 2-10 keV luminosities compared with the XRB model predictions 
(Gilli et al. 2007). Data points refer to this work (filled squares at low red- 
shift with upward-pointing arrows) and literature work (all the other sym- 
bols, at z 1 —2.5; reference papers are listed in the bottom part of the 
figure), and have to be compared with the corresponding color curves. 

since it is available for all of the 25 sources under consideration; 
however, any correction to the [O III] bias would decrease our esti- 
mate by ~ 25 per cent at most. 

Recently, R08 published an extended list of SDSS Type 2 
AGN candidates. They also reported new values for the [O III] line 
luminosity, which is measured using both a fitting with a Gaus- 
sian feature and with a non-parametric method (see Section 2.4 
of R08 for a detailed description), besides adopting the most up- 
to-date spectro-photometric calibration algorithm from the SDSS. 
Since our project was originally intended to provide an X-ray char- 
acterization of all of the sources of the Z03 sample in the quasar 
locus (i.e., at high X-ray luminosities, predicted accordingly to the 
[O III] vs. hard X-ray correlation of M94; see Fig. [Hi, in the fol- 
lowing we will keep the original [O III] line luminosity reported 
by Z03. However, it is interesting to analyze thoroughly what the 
main differences would be in case of source selection on the ba- 
sis of the [O III] line luminosity as reported in R08. If we adopt 
the same selection of logLjoni|> 9.28 L(-., using the R08 source 
catalog instead of the original catalog of Z03, we obtain a sample 
comprising 26 sources, 22 in common with the source list reported 
in Table[5] On the one hand, sources #16, #30 and #244 (two clas- 
sified as Compton thick on the basis of [O III] only, and one as 
intermediate/Compton-thin) would not be present anymore for our 
estimates of the volume density of SDSS Compton-thick quasars. 
On the other hand, sources #117, #141, #189 and #213 would be 
included. While source #117 must be removed because of its ra- 
dio loudness (its exclusion was described at the beginning of this 
Section), a similar destiny would be reserved to source #189, for 
the same reason. The other two sources lack any X-ray constrain 
to date, so they would not be useful for the following analysis. In 
the light of the limited number of changes implied by the adoption 
of the R08 catalog (and the way we account for the Compton-thick 



AGN predicted using only [O III], as described above), we are con- 
fident that the main conclusions of the present work (reported in 
the following paragraphs) will not change significantly. We also 
note that, for the sources above the chosen [O III] line luminosity 
threshold in both catalogs (22 objects), the logarithms of the lumi- 
nosities derived for [O III] differ, on average, by ~ 0.09, with no 
obvious evidence for systematic effects. 

In order to get the space density of Type 2 quasars above 
the chosen luminosity threshold, we integrated the SDSS Type 2 
quasars luminosity function published by R08 from logLjo iil| > 
9.28 Lq to the maximum luminosity for which they have data, 
i.e. logLjo ui]=10.3 Lq and logLjg ni] = 10.0 L;, for objects in the 
range z = 0.3 — 0.5 and z = 0.5 — 0.8, respectively. For the ob- 
jects in the z = 0.5 — 0.8 bin, we considered both curves in Fig. 5 
of R08, in which the measured space densities depend on a dif- 
ferent correction applied to the Type 2 quasar selection function 
(see R08 for details). These integrations return a space density 
of 1.2 X 10"** Mpc"^ for Type 2 quasars at z = 0.3 - 0.5 and 
( 1 .9 — 2.6) X 10"** Mpc"-' (depending on the correction mentioned 
above) for Type 2 quasars at z = 0.5 — 0.8. Since 5 out of 6 ob- 
jects in our z = 0.3 — 0.5 sample and 12 out of 19 objects in 
our z = 0.5 — 0.8 sample are Compton-thick candidates, we esti- 
mate the space densities of [O IIl]-selected SDSS Compton-thick 
quasars to be 1.2 x 10"^ x (5/6) = 10"^ Mpc"^ and (1.9-2.6) x 
10"^ X (12/19) = (1.2- 1.6) x 10"'* Mpc-^ at z = 0.3 -0.5 and 
z = 0.5 — 0.8, respectively. Based on the observed [O III] vs. intrin- 
sic L2-.\i^keV relations published by M94, the above space densities 
should then refer to Type 2 quasars with logLx > 44.6 erg s^'. In 
Fig.|4]we show the comparison between the space densities derived 
for our sample and those published for other samples of Compton- 
thick AGN candidates at different redshifts and luminosities. The 
densities presented in this work cover a corner in the luminosity 
vs. redshift plane of Compton-thick AGN so far unexplored. We 
also show the expected redshift evolution of the space density of 
Compton-thick AGN with different intrinsic luminosities based on 
the XRB model of Gilli et al. (2007). Our measurements are about 
a factor of 2-3 lower than the expectations of the XRB model, in 
which Type 2 Compton-thick quasars are assumed to be as abun- 
dant as unobscured quasars with the same intrinsic luminosity. A 
number of Compton-thick AGN smaller than expected (by a fac- 
tor of 2, corresponding to ~ 2a) might also be suggested by recent 
Integral (Beckmann et al. 2009) and Swift (Tueller et al. 2008) sur- 
veys of very local objects, which would point towards a revision 
of the XRB models (although this issue is still matter of debate, 
see Malizia et al. 2009). For instance, a lower space density of 
Compton-thick AGN features the recent model by Treister, Urry 
& Virani (2009). 

A few caveats have nonetheless to be kept in mind: indeed, 
R08 consider their estimates to be lower limits, since, if the NLR 
is extincted, a significant fraction of obscured AGN may escape 
detection based on narrow optical emission lines (other arguments 
are discussed in R08). In addition, the intrinsic X-ray luminosity of 
these Compton-thick Type 2 quasar candidates has been obtained 
using an observed [O III] vs. intrinsic X-ray luminosity relation de- 
rived by M94 for a local sample of Seyferts, which generally cover 
a lower [O III] luminosity range, and have not been selected on the 
basis of their [O III] flux or luminosity. As mentioned in ^2.31 other 
relations between [O III] and X-ray luminosities have been pub- 
lished by H05 and P06, but none of them is based on a purely [O III] 
selected sample, although H05 made an attempt to build a sample 
of [O IIl]-bright local Seyferts as a surrogate. For their sample of 
Seyfert 1 objects, i.e. for those sources in which the intrinsic X-ray 
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emission can be estimated directly, H05 derived X-ray luminosities 
on average a factor of 1.5 lower than those derived by M94 at a 
given [O III] luminosity [log Lx /L^oin] = l -59 vs 1.76]. Therefore, 
using the relation by H05 we would estimate for our sources intrin- 
sic luminosities of logLx > 44.4 erg s^'. In the considered XRB 
synthesis model, the space density of obscured AGN, and therefore 
also of Compton-thick AGN, is expected to increase with decreas- 
ing intrinsic luminosity. At z = 0.3 — 0.8, the XRB space density 
of Compton-thick AGN with logLx > 44 A erg s^^ is expected to 
be about an order of magnitude larger than that measured in this 
work. On the other hand, when using the relation by P06, which 
is essentially based on optically selected AGN samples, the [O III] 
luminosity threshold translates into an X-ray luminosity threshold 
of logLx > 44.8 erg s^'. At these luminosities, the model space 
density of Compton-thick Type 2 quasars is very close to our mea- 
surements. Summarizing, given the uncertainties mentioned above 
(see also § 12.31 ), at z = 0.3 — 0.8, Compton-thick Type 2 quasars 
selected by the SDSS are compatible with being as abundant as un- 
obscured quasars of similar luminosities (see also Alexander et al. 
2008). 



5 SUMMARY 

We have presented the most up-to-date results regarding the prop- 
erties of Type 2 quasar candidates selected from the SDSS and ob- 
served by Chandra. In particular, the objects presented in this work 
consists on two samples, one comprising 12 Type 2 quasars se- 
lected in Chandra Cycle 8 among the most luminous [O III] emit- 
ters, and the other comprising three additional candidates at lower 
[O III] luminosity retrieved from the Chandra archive. The full 
sample spans the redshift range z=0.40-0.73. The main results are 
summarized below: 

• Our Chandra exploratory observations (~ 10 ks for each point- 
ing) support previous findings that SDSS Type 2 quasars are X-ray 
faint: two sources were not detected by Chandra and nine of the de- 
tected sources have less than 10 counts in the observed 0.5-8 keV 
band. For only three sources it was possible to perform a low-to- 
moderate quality spectral analysis, providing evidence for column 
densities in the range Nu_ ^ lO-'^-lO^-' cm in the source rest 
frame. 

• The assumption of the [O III] luminosity as a proxy of the nu- 
clear emission, coupled with literature [O IIl]-hard X-ray correla- 
tions (with all the caveats of this approach being extensively dis- 
cussed in 32.31 and indicates that a significant fraction (about 
two-third) of the SDSS Type 2 quasars with logL|o ni|> 9.28 L0 
presented here and in V06 may be Compton thick (see also Ptak et 
al. 2006 and Lamastra et al. 2009). 

• Further indications of heavy obscuration in a sizable fraction 
of SDSS Type 2 quasars come from the analysis of archival Spitzer 
observations, which have been used to estimate the intrinsic (i.e., 
de-absorbed) hard X-ray AGN strength assuming some recent cor- 
relations between either the rest-frame 5.8iim or the 12.3(i.m lumi- 
nosity and the 2-10 keV luminosity. Using the combined [O III], 
mid-IR and X-ray information (by means of luminosity ratios; 
see Table |5]( it is possible to distinguish the Compton-thick from 
the Compton-thin Type 2 quasar population. Overall, we find that 
~ 50 per cent of SDSS Type 2 quasars with logLp ui]> 9.28 Lq 
appear to be obscured by Compton-thick material based on both 
the ix.meas^^x.mid-lR (where mid-IR corresponds to rest-frame 
12.3jim) and Lx.mea.s/i-x.[Olll| ratios. Moderate-quality X-ray spec- 
tra showing a strong iron Ka emission line over a flat continuum 



would provide a further, clean indication of Compton-thick absorp- 
tion. However, given the low X-ray fluxes of these sources, such 
spectral approach, extended to sizable samples, probably demands 
next-generation of sensitive X-ray detectors. 

• By integrating the SDSS Type 2 quasar luminosity function 
of R08 above logLjoni] > 9.28 Lq (corresponding to logLx > 
44.6 erg s^' under the assumption of the M94 relation; see Fig.lTJ, 
we obtain an estimate of the space density of Type 2 quasars above 
this luminosity threshold of 10"** Mpc"-* and (1.2 — 1.6) x 
10"^ Mpc"3 at z = 0.3 -0.5 and z = 0.5 -0.8, respectively. At face 
value, these space densities are a factor 2-3 lower than the expec- 
tations of the XRB model by Gilli et al. (2007) and would suggest 
that Compton-thick Type 2 quasars are less abundant than unob- 
scured quasars (Treister et al. 2009). However, given the uncertain- 
ties in the SDSS Type 2 quasar luminosity function (mainly related 
to possible extinction within the NLR) and in the assumption of the 
[O IIl]-hard X-ray correlation (see we can conclude that the 
estimated space density of SDSS Compton-thick Type 2 quasars 
at z = 0.3 — 0.8 is likely a lower limit. In other words, in this still 
poorly explored redshift range (recent results on space densities of 
heavily obscured quasars are mostly at z ~ 1 — 2.5), Compton-thick 
Type 2 quasars are consistent with being as abundant as unobscured 
quasars of similar intrinsic luminosities. 
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